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Abstract. The Coma cluster of galaxies was observed with XMM-Newton in 12 partially overlapping pointings. 
We present here the resulting X-ray map in different energy bands and discuss the large scale structure of this 
cluster. Many point sources were found throughout the observed area, at least 11 of them are coincident with 
bright galaxies. We also give a hardness ratio map at the so far highest angular resolution obtained for a cluster 
of galaxies. In this map we found soft regions at the position of bright galaxies, little variation in the central 15 
arcmin, but some harder regions north of the line NGC 4874 - NGC 4889. 
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1. Introduction 

The Coma cluster of galaxies was one of the performance 
verification targets of XMM-Newton. The aim was to ver- 
ify XMM-Newton's ability to map large extended X-ray 
sources. This is not an easy task because errors can occur 
in several areas. 

The telescope vignetting plays a major role when over- 
lapping the same region of an extended source, measured 
at different off-axis angles. Because of the vignetting, the 
surface brightness of a region measured off axis seems to 
be lower than the same region measured on axis. The 
reason for that is the increasing loss of mirror collecting 
area with increasing off-axis angle, which in addition is 
energy dependent. By de- vignetting the different obser- 
vations and normalizing the observations to the obser- 
ving times the surface brightnesses of the same source 
regions are corrected and must show the same number 
of counts/arcmin^/sec. If done correctly, detector-edges 
or regions of detector gaps should not be visible as inho- 
mogenities in the surface brightness of the merged obser- 
vations. 
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The backgrounds in the image are problematic. The 
part of the background that is produced by particles inside 
the detector is not distinguishable from X-rays, is not vi- 
gnetted, and is in general distributed homogeneously over 
the detector face (except for the lower end of the spec- 
tral range and around the detector intrinsic emission lines 
around 8 keV, see Briel et al. 200C| ). Hence this compo- 
nent must be subtracted before de-vignetting is applied. 
The extragalactic X-ray background is vignetted in the 
same manner as diffuse emission from an extended source 
and may therefore be de-vignetted in the same manner. 
But there is a third component in the background, the 
low energy protons, that are again indistinguishable from 
X-rays. They come through the mirror system and show a 
different vignetting compared to X-rays (see Striider et al. 
pool ). In addition, their spectrum and their fiux is time 
variable. Because of the different vignetting, it is essential 
to subtract a "proton image" from the measured image 
before de-vignetting is applied, especially when the sur- 
face brightness of the proton induced background is in 
the same order of magnitude as the surface brightness of 
the extended source. 

Coma was chosen because it is a well studied source 
at all wavelength ranges. With a diameter of the main 
cluster, which is several times larger than the FOV of the 
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Table 1. Journal of Observations 



Name of pn-Camera observing times (ksec) 

Date Observation RA(2000) DEC(2000) planned performed effective 
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EPIC cameras of about 26 arcmin, it is necessary to mo- 
saic the cluster with XMM-Newton if the large scale struc- 
ture and dynamics of the cluster are to be investigated. 
Although Coma was long believed to be the archetype of 
a relaxed cluster of galaxies, the X-ray image taken dur- 
ing the Rosat All Sky Survey revealed the complex large 
scale structure, especially the subgroup around NGC 4839 
probably falling into the main cluster (Briel et al. |l992|) . 
Subsequent pointed observations with Rosat showed that 
Coma has irregular structure on different scales and was 
most probably formed by hierarchical clustering of several 
subunits, whi ch ar e still visible in their X-ray emission 
(White et al. 1993 ). Observations with ASCA, and also 
the pointed Rosat observation, showed evidence of tem- 
perature structure in Coma (Honda et al. 1996; Briel & 
Henry |1998| ; Donnelly et al. |1999| ; Watanabe et al. |1999D . 
The interpretation of the temperature structure favors the 
ongoing merging of the subgroup around NGC 4839 with 
the main cluster, a scenario which is still debated. 

In this letter we will present the result of merging 
12 partially overlapping observations of the Coma cluster 
with the EPIC pn-CCD-detector on board XMM-Newton. 
In Section 2 we will describe the observations and present 
the data reduction and the method of merging the indi- 
vidual observations. In Section 3 we will discuss some of 
the many point sources in the image, the large scale struc- 
ture of the X-ray surface brightness distribution, and the 
hardness ratio map of the inner cluster region. 



2. Observations and Data Reduction 

The region of the Coma cluster of galaxies was observed 
by the XMM-Newton Satellite (Jansen et al. ^00l[ ) during 
the performance verification phase. In Table 1 we show 
the journal of the observations with information on the 
pointing directions and observing times. 



be well-calibrated for spectral analysis. However due to 
operational difficulties, less data coverage of the Coma re- 
gion was obtained in these cameras. The pn-CCD camera 



(Striider et al. 2001) was operated in the extended full 
frame mode (see the above papers for a description of the 
different modes). This latter mode provides potentially the 
best imaging performance of the pn-CCD camera, but at 
this date its spectral calibration is less reliable than that 
of the MOS observation. The analysis of th is pap er there- 
fore contrasts with that of Arnaud et al. ( 2001 ), in that 
we concentrate more on the gross imaging properties and 
less on the spectral imaging. The hardness ratio map that 
we present, based on flux estimate in wide energy bands, 
is not sensitive to the remaining PN calibration uncertain- 
ties. For all observations, the medium filter was used (see 
2Ci0l| ). 



Turner et al. 



The two MOS-CCD cameras (Turner et al. |200l|) were 
operated in their full-frame mode, which we believe to 



During the observations the proton induced back- 
ground often varied, and in flares reached count rates up 
to several hundred counts/sec in the energy band from 
0.15 to 10 keV. To increase the signal-to-noise ratio and 
also to minimize the above described problems of the de- 
vignetting of residual proton induced events, we binned 
the 10-15 keV events in 100 sec time bins. From this light 
curve, times of bins with count rates above 100 counts/100 
sec were identified and excluded from further processing. 
We have chosen the high energy region because we do not 
expect many events from Coma in that band. This rate 
criterion gave in all but Coma 4, 6, and 7 reasonably long 
effective exposure times, which are given in Table 1. The 
rate criterion for Coma 4 and 7 were 200 counts and 800 
counts respectively. We did include these high background 
times in order to get the mosaic of Coma as complete as 
possible. We do have to be careful with the interpretation 
of those two fields. Since the Coma 6 field covers part of 
the center of the cluster, we did apply the stringent rate 
criterion, therefore 3.7 ksec exposure time is usable, which 
is only 25% of the total time pointed at this position. 
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Fig. 1. The merged EPIC-pn image of the Coma cluster of galaxies in the 0.3 to 2.0 keV energy band. Indicated are 
the tentative identifications of point sources. 



In addition to the above screening, we also excluded 
"warm" pixels from the images by checking each indivi- 
dual observation in detector coordinates and searching for 
single pixels with statistically significant higher numbers 
of events compared to adjacent pixels. 

Images were produced in three different energy bands: 
0.3 - 2 keV, 2 - 5 keV and 5 - 7.15 keV. We stopped at 7.15 
keV in order to avoid the strong detector intrinsic emis- 
sion lines of Ni-K, Cu-K and Zn-K, spread over approxi- 
mately 2 keV. At a later time, when we learn more about 
the background, the higher energies will be included. An 
other handicap was the fact that, at the time of analysing 
the pn-data the correction files for the charge transfer inef- 
ficiency (CTI) were only available for the full frame mode 
(not the extended full frame mode used for our obser- 
vations), so they were used to correct for the CTI. This 
effects the corrected amplitude of the events at the level of 
5-25 eV depending on the energy and the position in the 
detector. But for producing images in the above bands, 
this systematic error is negligible. All valid events were 
used in the images. Valid events are all single events plus 
recombined split events, whose position is defined as the 
position of the pixel containing the maximum amplitude 



(see Dennerl et. al. 1999 for a detailed explanation). To 
minimize the gaps between CCDs, events at the border 
pixels are also included, although their amplitude could be 
slightly under estimated because of a possibly not recorded 
split companion. Again, for making images in wide energy 
bands, this systematic error is negligible. 

In the same energy bands, and with the same se- 
lection criteria, background images were produces from 
a dedicated "background event file". This file contains 
132.8 ksec of pn-data accumulated from 5 other observa- 
tions, like the Lockman Hole (Hasinger et al. |200l| ), con- 
taining no extended X-ray sources and only weak point 
sources. The single observations are between 13 and 33 
ksec long, spread over the 3 months CAL/PV phase of 
XMM-Newton. These data were also screened for flares 
and warm pixels, and obvious point sources were taken 
out. Comparing the single observations of the screened 
event files, the background rate in the energy band from 
0.15 to 10 keV is very stable, the systematic difference be- 
tween the minimum and maximum rate is 8% of the mean. 
From each Coma image the corresponding background im- 
age, normalized using the effective observing times, was 
subtracted. Even though we restrict ourselves to observing 
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Fig. 2. The merged EPIC-pn image of the Coma cluster of galaxies in the 2.0 to 5.0 keV energy band. 



times with "low" background, at a distance of 30 arcmin 
from the center of Coma the background is compareable 
to the surface brightness of Coma. At the center of Coma 
the background is ^ 1.5% and 2.5% of Coma's surface 
brightness in the bands 0.3 to 2 keV and 2 to 7.15 keV 
respectively. 

In a further step, exposure maps were produced for 
the different energy bands, also normalized to the corre- 
sponding effective observing times. These exposure maps 
contain the energy dependent vignetting. 

The individual images in the lowest energy band were 
searched for point sources, and the positions of point 
sources in overlapping regions were then used to ahgn and 
merge both the background-free X-ray images and the ex- 
posure maps. In a last step, the merged background free 
X-ray images were divided by the corresponding merged 
exposure maps. The result is the surface brightness distri- 
bution of the Coma field in the different energy bands. In 
the Figures 1 and 2 we show the 0.3-2.0 keV image and 
the 2-5 keV image of Coma taken with the pn camera. 

To visualize better the low surface brightness of the 
cluster we smoothed the images with two-dimensional 
Gaussian filters of variable sizes (with a sigma from 9 
arcsec to 96 arcsec). We used these images to look for 



spectral variations in the central Coma region. A hard- 
ness ratio map was produced by calculating the im- 
age HR = (image(2-7.15keV)-image(0.3-2keV))/(image(2- 
7.15keV)-|-image(0.3-2keV)). In Fig. 3 we show the inner 
25.5x23.3 arcmin^ of this hardness ratio map. The over- 
lain contours are from the smoothed surface brightness 
distribution in the 0.3-2 keV band. We performed simu- 
lations with XSPEC, using a Raymond & Smith model 
(with the fixed parameters z — 0.0232, abundance = 0.25 
and Nh = 8.95 10^^ cm-^ from Dickey & Lockmann [l990| ) 
and taking into account our current understanding of the 
errors in the preliminary calibration of the extended full 
frame mode. With these simulations we converted the 
hardness ratios into " temperatures" . The resulting con- 
version is indicated as color bar with the corresponding 
temperatures in Fig. 3. 

The dynamic range of the hardness ratio in the map 
is of the order of 20%, corresponding to a "temperature" 
variation of ~ 5 keV. A systematic error of the background 
model of 10% would be neglectible in the inner 15 arcmin 
diameter. At larger radii the ratio would be different by ^ 
3%, corresponding to a "temperature" difference of about 
1 keV. 
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Fig. 3. The hardness ratio map of the inner Coma region in the energy bands from 0.3-2 keV and 2-7.15 keV. The 
overlain contours are from the smoothed surface brightness distribution in the energy band from 0.3-2 keV 



3. Discussion 

The Coma cluster is the nearest very rich cluster of galax- 
ies. It is probably the best studied cluster at all wave- 
lengths and is one of the brightest extragalactic X-ray 
sources in the sky. Until now, the deepest X-ray obser- 
vation of Coma is that of White et al. ( 1993| ) consisting of 
4 exposures of about 19 ksec each by the ROSAT PSPC. 
The EPIC pn observations reported here are of compara- 
ble length, thus are ^8 times deeper, while simultaneously 
extending over '-~^3 times wider energy range and with ~4 
times better angular resolution. 

Perhaps not surprisingly given the better spatial re- 
solution and deeper observation, the most noticable dif- 
ference between our XMM images in Figures 1 and 2 
compared to Figure 2 of White et al. (1993) is the myr- 
iad of newly detected point-like sources. Dow & White 
( |1995D discuss the ^^25 point sources in the ROSAT ob- 
servations of Coma. There are at least 75 point sources in 
the smaller solid angle we have analyzed. We have made 
a first attempt to identify these sources by overlaying our 
X-ray map on an optical image of the region and then 
consulting NED when there appears to be a coincidence. 



Many bright galaxies seem to be X-ray sources including: 
NGC 4827, 4839, 4858, 4860, 4874, 4889, 4911, and 4921, 
and IC 4040 and 4051. However some faint galaxies, with 
magnitudes from 18 to 20, could also be identifications 
based on positional coincidence. At least three quasars, 
QSO 1256-^281, QSO 1258-^280, and QSO 1259+281, are 
detected. However, the bulk of the sources are unidentified 
in our initial analysis. 

The NGC 4911 and NGC 4839 groups contain concen- 
trations of point sources at our s ensitivity. Dow & White 
(1995) and Neumann et al. ( 2001 ) discuss the particularly 
interesting morphology of the NGC 4839 group. 

The large-scale diffuse emission in the soft band shown 
in Figure 1 is very similar to the ROSAT observations of 
White et al. ( |1993| ) in the same band. The X-ray emission 
of the Coma cluster is lumpy near its center but becomes 
smoother at larger radii. An additional lump, apparent 
in the ROSAT observation but not discussed by White 
et al. (|1993|) , is directly NE of NGC 4839. It may be yet 
another group on its way into the main cluster to the ENE. 
This lump may be preceeding the NGC 4839 group from 
the direction of A1367 since calculations of cluster growth 
indicate that the groups that merge to form clusters flow 



6 



U.G.Briel et al.: A Mosaic of the Coma Cluster of Galaxies with XMM-Newton 



along filaments connecting the clusters. The large-scale 
diffuse hard emission is very similar to the soft, as Figure 
2 shows. 

Because the calibration of the spectral response of the 
extended full frame mode of the pn-CCD is not yet avail- 
able, we do not yet want to perform spatially resolved 
spectral fits. We can, however, make a hardness ratio map 
using wide bands, which we show in Figure 3. While these 
maps sacrifice spectral details, such as temperatures and 
abundances, they are the highest angular (~40 arcsec) and 
spatial (^26 kpc for Hq = 50 km Mpc~^) resolution 
spectral maps ever obtained for a cluster of galaxies. 

Using the conversion from hardness ratio to tempera- 
ture, given in Fig. 3, we have checked our resulting mean 
temperature in some regions which are known from other 
measurements. For the over-all temperature within a ra- 
dius of 10 arcmin, centered on NGC 4874, we determined 
a temperature of 8.8 ke V. Th e GINGA value is 8.21 ± 
0.09 keV (Hughes et al. |1993D , the ASCA value is 9. ± 
0.6 keV (Donnely et al |1999[ ) a nd the XMM-MOS value is 
8.25 ± 0.10 keV (Arnaud et al |200l[ ). We agree to better 
than 10% with these previously determined values, which 
would be expected given our understanding of the uncer- 
tainties of the calibration and of the background model. 
In a further check, we determined the temperature around 
NGC 4874 within a radius of 0.5 arcmin and within a ring 
from 1 to 5 arcmin. Our result is 7 and 9.2 keV, the MOS 
result (Arnaud et al |200l|. Fig. 8) is 6.6 and 8.6 keV re- 



spectively, again in agreement within 10%. We found the 
same kind of agreement on other places, taking into ac- 
count the large area of 3.5 x 3.5 arcmin^, used for the 
MOS data by Arnaud et al. ( |200lD . 

The large-scale descrip tion of the data in Figure 3 
agrees with Arnaud et al. (2001): there is cool gas to the 
south-east and hotter gas to the south-west and west. 
It agrees to so me ext ent with Donnely et al. ( |l999|) , 
Watanabe et al. ( |l999| ) and Briel & Henry (|l998|) inves- 
tigations. In particular the south-east cool region was al- 
ready identified by Donnely et al. ( |1999[ ) and a hot region 
in the south-west /west direction, but somewhat further 
away, was found by Watanabe et al. ( |l999[ ). 

There are spectral variations on the smallest scales 
that we can measure. Softer emission is associated with 
many compact sources some of which are associated with 
galaxies as discussed previously. Although the smoothing 
of the original surface brightness distribution smears out 
the point sources to some extend, a significant drop in the 
radial temperature profile within 1 arcmin of NGC 4874 
was also seen by Arnaud et al. (2001), Fig. 8. The Rosat 
temperature map of Briel & Henry (1998) also shows 
cooler regions around the two central galaxies (although 
other cool regions in the center of their map are not con- 
firmed). If these soft regions, associated with galaxies, are 
extended, we can postulate three possible explanations. 
First, they are local potential wells centered on cluster 
galaxies as is was discussed by Vikhlinin et al. (1994). 
The shallower potential wells could only trap lower tem- 
perature gas. Second, they are the interstellar emission 



that the cluster galaxy has been able to retain. Third, 
integrated emission from low mass X-ray binaries which 
Irwin & Sarazin (1998) have discussed as the source of soft 
emission from elliptical galaxies. 

Perhaps the most interesting small scale spectral struc- 
tures are the isolated hard spots that do not appear to be 
associated with any obvious point-like X-ray emission. For 
the most part they are not associated with an optical ob- 
ject in NED. The five hard spots north of the NGC 4889 
- NGC 4874 line may be the cause of the hot region in the 
same general location detected by Donnelly et al. (1999) 
in their ASCA observations. It seems that these hot spots 
are wa shed out in the temperature map of Arnaud et al. 
(2001) because of the larger (3.5x3.5 armin^) integration 
boxes. They might have been seen however in the even 
larger integration region of ASCA, if they are associated 
with variable sources like QSO's. But we can not entirely 
eliminate an instrumental cause for some of them because 
some are aligned parallel to the CCD orientation and on 
the boarder of the chips. 
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